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SUMMARY 

An experimental study of low Reynolds number nozzle flow was performed. 

A brief comparison was made between some of the experimental performance data 
and performance predicted by a viscous flow code. The performance of 15°, 20°, 
and 25° conical nozzles, bell nozzles, and trumpet nozzles was evaluated with 
unheated nitrogen and hydrogen. The numerical analysis was applied to the con- 
ical nozzles only, using an existing viscous flow code that was based on a 
slender-channel approximation. Although the trumpet and 25° conical nozzles 
had slightly better performance at lower Reynolds numbers, It Is unclear which 
nozzle Is superior as all fell within the experimental error band. The numer- 
ical results were found to agree with experimental results for nitrogen and for 
some of the hydrogen data. Some code modification Is recommended to Improve 
confidence In the performance prediction. 


INTRODUCTION 

Electrothermal thrusters, such as arcjets and reslstojets, have the poten- 
tial for high performance. These thrusters operate at high temperature and low 
thrust levels, which result In low Reynolds number nozzle flow. In the low 
Reynolds number regime, the flow boundary layer occupies a large portion of the 
flow area. These thick boundary layers can result In reduced thruster perform- 
ance because of Increased viscous losses. It Is felt that through proper 
selection of the nozzle geometry these losses could be minimized. 

In the past, work has been done In the low Reynolds number regime for 
supersonic nozzles (refs. 1 and 2), for chemical laser nozzles (refs. 3 to 5), 
and for low-density hypersonic nozzles (refs. 6 and 7). General consensus from 
these studies was that the viscous boundary layers In the nozzles were very 
large and that the flow could become fully viscous, with no Invlscld core, as 
the Reynolds number was decreased. Rothe (ref. 1), In his electron-beam 
studies, showed that In supersonic nozzles, at a Reynolds number of 500, a 
small Invlscld core exists In the flow and that below a Reynolds number of 300 
the flow was fully viscous with no Indication of an Invlscld core. Driscoll 
(ref. 3) conducted a boundary-layer study In chemical laser nozzles and found 
that for axlsymmetrlc nozzles the flow becomes fully viscous below a Reynolds 
number of about 1000. Whitfield (ref. 6) performed an Investigation Involving 
hypersonic nozzles and found that the boundary layer occupied as much as 
80 percent of the nozzle exit area. These thick boundary layers contributed 
to the viscous losses In the nozzles by displacing the Invlscld core, so that 
In the divergent section, the core density was Increased and the velocity was 
decreased. 

Early work In the area of nozzle performance has generally dealt with con- 
ical nozzles of varying half angles and area ratios. Murch, et al . (ref. 8), 


examined 10°, 20°, and 30° conical nozzles and also tested bell and horn noz- 
zles. The tests were performed with heated hydrogen and nitrogen for Reynolds 
numbers ranging from 400 to 4000. The results showed performance losses of up 
to 30 percent at low Reynolds numbers. A more recent study by Brophy, et al. 
(ref. 9), compared a modified bell and a 19° conical nozzle for arcjets. 
Recently, Grlsnlk, et al . (ref. 10), Investigated the flow through four dif- 
ferent nozzle contours. This study was performed with unheated hydrogen and 
nitrogen, and Indicated that within experimental error the four contours had 
essentially the same performance. 

Studies by Edwards and Jansson (refs. 11 to 13) Investigated the loss 
mechanisms of a reslstojet. Boundary-layer effects on nozzle performance and 
the Importance of operating Reynolds number In the design of reslstojet nozzles 
were discussed. Concerns have also been raised about difficulties In measuring 
the performance of thrusters In vacuum environments. Yoshlda, et al. (ref. 14), 
showed that the ambient pressure of the vacuum chamber has a significant effect 
on performance measurements. Recently, Manzella, et al. (ref. 15) found that 
the effects on performance measurements occurred when testing with hot gases 
and suggested that these effects were due to convective heat transfer losses. 

The study presented herein used unheated gas flow to eliminate this loss 
mechanism. Another study by Sovey, et al. (ref. 16), showed that the ratio 
between the nozzle Inlet pressure and the ambient pressure In the vacuum 
chamber needed to be lower than 10~ 3 to assure that a shock was not present 
In the nozzle and thus affecting performance of the nozzle. This consideration 
becomes particularly Important for measurements at higher flowrates. 

Several studies have examined the transonic region of supersonic nozzles. 
Back, et al. (refs. 17 and 18), Investigated the Influence of the contraction 
section of conical nozzles. Moreover, they measured and predicted flows. The 
performance was found to be relatively Insensitive to the nature of the flow 
In the convergent section, but the radius of curvature through the throat 
region did affect the discharge coefficient. Work by Campbell and Farley 
(ref. 19) on converging-diverging nozzles Indicated that a decrease In the 
divergent half angle resulted In Increased thrust. Further, Back and Cuffel 
(ref. 20), found that for a relatively high Reynolds number flow, (10 6 ), the 
mass flow through the throat was affected by the throat configuration. The 
discharge coefficient was shown to decrease for a decrease In the ratio of 
throat radius of curvature to throat radius. This was corroborated by Cuffel, 
et al. (ref. 21), and again the convergent half angle was found to have little 
Influence on discharge coefficient. Hopkins and Hill (ref. 22) developed a 
method to numerically predict the flowfleld In the transonic region of nozzles. 
It was shown that the most significant geometric factor Influencing flow 
through the throat was the radius of curvature. Apparently, the shape and 
location of the sonic line was altered by the curvature of the throat, and the 
effect of the convergent angle was found to be secondary, unless the radius of 
curvature was less than the throat radius. Kuluva and Hosack (ref. 23) exam- 
ined the discharge coefficient In low Reynolds number flow. Their work Indi- 
cated that velocity slip at the wall was significant at Reynolds numbers below 
10 3 and that the throat curvature became Important at Reynolds numbers below 
200. Milligan (ref. 24), also Indicated that slip velocity at the nozzle wall 
should be taken Into account. 

Rae (refs. 2 and 25) developed a numerical code which used a slender- 
channel approximation to predict viscous effects and the performance of 
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low-density nozzle flows. Rae's method used equations for the slender-channel 
approximation developed by Williams (ref. 26). The latter was an analytical 
study which Involved compressible flows In convergent-divergent channels and 
yielded results regarding viscous effects, pressure ratio at the throat and 
discharge coefficient. Rae used the slender-channel approximation In conjunc- 
tion with slip boundary conditions to predict nozzle performance In most of the 
Reynolds number range appropriate for electrothermal thrust devices. However, 
this technique was only capable of obtaining solutions for certain combinations 
of Reynolds number, wall angle, wall temperatures, and gas properties. The 
envelope of conditions for the successful use of this code was not fully 
defined. Results obtained from the code have been compared to Navler-Stokes 
solutions by Mltra and Fleblg (ref. 27). These solutions showed differences 
originating at the throat, but the solutions converged farther downstream. 

The work presented here Includes an experimental study of low Reynolds 
number nozzles and a brief comparison of some of the experimental performance 
data to numerical performance predictions using the slender-channel approxima- 
tion. In attempting to Improve the performance of electrothermal thrusters, 
the evaluation of the performance of various nozzle contours was considered to 
be significant. Nozzle performance Is affected by many factors Including noz- 
zle contour, area ratio, and operating temperature. In the experimental por- 
tion of this Investigation, several different nozzle contours were selected for 
performance evaluation with unheated gas flow. It was felt that loss mecha- 
nisms could be better quantified and Identified using unheated gases. 

The series of experiments compared the performance of conical nozzles with 
15°, 20°, and 25° half angles, bell nozzles and trumpet-shaped nozzles. Oata 
with unheated nitrogen and hydrogen were taken for each nozzle. The nitrogen 
data Included Reynolds numbers up to 6000 and hydrogen data were for Reynolds 
numbers up to 3500. The variation of specific Impulse efficiency and discharge 
coefficient for each nozzle as a function of Reynolds number was evaluated for 
different area ratios and nozzle contours. 

In addition to the experimental study, a numerical analysis of the conical 
nozzles was undertaken. This analysis used an existing code developed by Rae, 
discussed previously. The code was run In an attempt to determine whether the 
slender-channel approximation was valid for conical nozzle designs and flow 
conditions considered here. 


NOMENCLATURE 

A CO de flowrate generated by the viscous flow code 
Agiven flowrate based on experimental data 
B Reynolds numbers based on reservoir conditions 

Cf thrust coefficient 

M* throat Mach number at centerline of nozzle 

M e nozzle exit Mach number at centerline of nozzle 
Re Reynolds number based on throat conditions 
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r/r* ratio of throat radius of curvature to throat radius 
X max axial location of nozzle exit normalized to throat radius 

X 0 axial location of convergent section normalized to throat radius 

c area ratio 

01 convergent half angle of nozzle 

02 divergent half angle of nozzle 


APPARATUS AND PROCEDURE 
Experimental Apparatus 

The facility used for the nozzle testing was a 0.91 m (3.0 ft) diameter 
by 1.82 m (6.0 ft) long vacuum tank. A 30 cm pipe connected the tank to a 
blower with a rated performance of 1x10 s llter/mln (3900 ft 3 /m1n (cfm)), run In 
series with a roughing pump with a rated performance of 2xl0 4 llter/mln 
(760 cfm). The facility was capable of pressures as low as m1d-10~* torr. 
Typical tank pressures with flowing gas varied from 6xl0 -3 to 0.2 torr, 
depending on gas type and mass flowrate. The tank pressure was measured using 
a cold cathode Ionization gauge, which had errors of up to ±20 percent In the 
low pressure range ( 1 0 -3 torr). Tank pressures less than 10 -2 torr were 
obtained In the tests with Reynolds numbers less than 2000. These tank or 
background pressure readings were used to correct all the measured thrust read- 
ings, and to determine the pressure ratio between the Inlet pressure and this 
background pressure. As mentioned above, a previous study (ref. 16) determined 
that this pressure ratio must be lower than 10 -3 torr In order to assure that 
a normal shock did not appear In the nozzle. 

The mass flowrate measurements used laminar flow, heated wire type flow- 
meters. For these flowmeters, the gas flow was proportional to wire Input cur- 
rent, for a constant wire temperature. Two flowmeters, covering the ranges 
from 0 to 5 standard llters/mln (slpm) and from 0 to 50 slpm, were used In 
these tests. The flowmeter was factory calibrated for hydrogen. The nitrogen 
values were obtained by using the manufacturer-supplied correction factor. The 
error In the flow measurement In the 0 to 5 slpm flowmeter was less than 
1 percent over the entire range. The flow measurement error of the 0 to 
50 slpm flowmeter was from 3 to 5 percent over the test range (6 to 26 slpm). 

The thrust measurements were obtained using a flexure-type thrust stand, 
capable of measuring forces as low as 2.2 mN (0.5 mlb). A schematic diagram 
of the thrust stand Is shown In figure 1. The thrust stand consisted of a hor- 
izontal mounting plate supported by two sets of four flexure plates. This 
arrangement allowed motion In the horizontal direction. The design of the 
flexures resulted In a large displacement for a small force. This design also 
reduced thermal and vibrational problems. The thermal expansion of the thrust 
stand components, when running hot flow tests, was negligibly small compared 
to the thrust stand displacement. Any high-frequency vibrations were prevented 
from reaching the thruster mounting plate due to the low resonant frequency of 
the flexures. A magnetic damper assembly was used to remove the low-frequency 
vibrations. The damper consisted of a permanent magnet attached to the 


4 



thruster mounting plate, and Inserted Into an annular electromagnet. The out- 
put of the electromagnet was dependent upon the movement of the thruster mount- 
ing plate. The movement of the thruster plate was detected by a linear 
variable displacement transducer (LVDT). The signal from the LVDT was sent to 
a chart recorder, which recorded the displacement. The LVDT output was also 
amplified and used to energize the electromagnet of the damper assembly. The 
magnetic field produced In this way then applied a force on the permanent mag- 
net opposing the motion of the thrust plate. 

The thrust range of the thrust stand could be varied by changing the 
stiffness of the flexures or the stiffness of the propellant feed tube, which 
was designed to also be a flexure. The feed line was 0.318 cm (0.125 In.) thin 
wall stainless steel tubing which was connected perpendicularly to the axis of 
motion of the thrust stand, and thus acted as an additional flexure. 

Calibration was performed by adding known weights to a pulley assembly and 
deflecting the thrust stand a given amount. The calibration was recorded on a 
chart recorder. The variability In thrust stand movement was small and repro- 
ducible. The thrust stand measurement variability was less than 5 percent at 
low thrust values and less than 2 percent at higher thrust levels, based on 
these calibrations. 

The effect of circulating gases In the vacuum chamber was determined by 
flowing gas through an orifice mounted close to the stand and checking for 
deflection of the thrust stand. The deflection ranged from 2.2 to 8.9 mN 
(0.5 to 2 mlb) depending on flowrate. The influence of circulating gases was 
negligible after Installing a wind shield between the nozzle exit and thrust 
stand. 

The nozzles, designed to be Interchangeable, were clamped to an assembly 
mounted on the thruster mounting plate. The assembly consisted of a 1.27 cm 
(0.5 In.) o.d. aluminum tube with stagnation pressure and temperature taps, 
feed line attachment and a flange which mated to the flange on each nozzle. 

The nozzles were clamped to this assembly as shown In figure 2. The gas inlet 
pressure was measured using a stainless steel, diaphragm tvpe transducer. The 
range of the pressure transducer was from 0 to 2.4xl0 5 N/m* (0 to 20 psla). 

The pressure tap was a 0.159 cm (0.0625 In.) tube located approximately 5.08 cm 
(2 In.) upstream of the nozzle throat. The pressure reading calibration was 
checked before each test and the error was found to be less than 1 percent of 
full scale. The gas temperature was measured using a stagnation type probe 
with a chromel-alumel thermocouple mounted In the gas stream about 5.08 cm 
(2 In.) upstream of the nozzle throat. Temperature measurement errors were 
estimated to be ±2 °C. 

Five sets of converging-diverging nozzles were designed for this study. 
Each set of nozzles had a different diverging contour with area ratios (exit 
area/throat area) of 25, 50, 100, 150, and 200. Each of the nozzles was 
machined from stainless steel and had a flange piece electron beam welded to 
It. The converging section was the same for each nozzle. It had an o.d. of 
1.27 cm (0.50 In.), an l.d. of 0.94 cm (0.37 In.), and converged at a 45° angle 
to a 0.152 cm (0.060 In.) throat diameter. 
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Three sets of conical nozzles were machined with 15°, 20®, and 25° half 
angles In the diverging section. Figure 3 provides an Indication of the range 
of sizes of the 20° conical nozzles. The 15° and 25° nozzles are similar, with 
slightly different diverging section length. A 15° conical nozzle with an area 
ratio of 25:1 was not available for testing. The design drawing of these noz- 
zles Is shown In figure 4, with the dimensions and area ratios. The throat 
section of these nozzles was a cylindrical section, 0.152 cm (0.060 In.) In 
length, with the radius of the converging and diverging sections machined so 
there was no corner or step. 

The fourth set of nozzles had a diverging section with a bell contour. 

The set of bell shaped nozzles Is shown In figure 5. These nozzles were 
designed using a technique developed by Rao (ref. 29). This technique assumed 
Isentroplc flow and used a variational method to design the wall contour. The 
method of characteristics was used to construct the nozzle contour for given 
flow conditions. The computer code used to develop the different bell contours 
required area ratio, radius of curvature Into and out of the throat, and the 
flow conditions. The result was coordinates that were used to plot the nozzle 
contour. The coordinates were subsequently supplied to a computer controlled 
lathe which machined the nozzles. Figure 6 shows the basic design of the bell 
nozzles as well as nozzle dimensions and area ratios. 

The final set of nozzles designed had trumpet shaped diverging contour 
sections. Figure 7 depicts the trumpet nozzle size variation with area ratio. 
This contour was designed as shown In figure 8. A circular arc was drawn 
starting from the throat and extending out to the diameter necessary for the 
different area ratios. These nozzles were also machined by using a computer 
controlled lathe. The dimensions and area ratios are listed In figure 8. 


Experimental Procedure 

Testing of the nozzles Involved consideration of four parameters: the 

nozzle contour, the nozzle area ratio, the gas type, and the Reynolds number. 

A total of 24 nozzles was tested. The nozzles with 15° half angle diverging 
sections had area ratios of 50:1, 100:1, 150:1, and 200:1. The 20° and 25° 
half angle conical nozzles, the bell contour nozzles and the trumpet shaped 
nozzles consisted of five nozzles each with area ratios of 25:1, 50:1, 100:1, 
150:1, and 200:1. Testing was done with two gases, hydrogen and nitrogen. The 
Reynolds number was based on throat conditions and was equal to four times the 
mass flowrate divided by pi, the viscosity and the throat diameter. The 
Reynolds number was varied by changing the mass flowrate. The range of 
Reynolds numbers tested was from 150 to 3500 for hydrogen, and from 500 to 6000 
for nitrogen. Two ranges of thrust were required to cover this range of 
Reynolds numbers: 3.0 to 25.3 mN (0.7 to 5.7 mlb) and 12.0 to 91.7 mN (2.7 to 

20.6 mlb). 

Prior to testing, the flexures on the thrust stand were adjusted to accom- 
modate the selected thrust range. The nozzle was then clamped In place on the 
thruster assembly (see fig. 2), the nozzle exit was plugged, and the assembly 
was pressure leak checked. If the gas pressure did not drop In approximately 
1 min, the system was considered leak tight. After the leak test, the thrust 
stand assembly was put Into the tank and the tank was sealed and evacuated. 

The gas feed system was also evacuated and purged several times to assure that 
the system was clean and only contained the test gas. When the tank pressure 
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stabilized, the thrust stand was calibrated. The calibration was performed at 
least twice prior to the recording of any thrust measurements to assure repro- 
ducibility, and to allow adjustment of the zero point If necessary. A mass 
flowrate was then set and the thrust zero checked before and after each thrust 
measurement. The thrust measurement was recorded twice for each flowrate (or 
Reynolds number). Data were taken for hydrogen flowrates from 1 to 26 slpm and 
for nitrogen flowrates from 0.5 to 6 slpm. This procedure was followed for all 
nozzles, with both gases and over both thrust ranges. In all tests, the gas 
was at ambient temperature. 

The following data were collected for each test point: flowrate. Inlet 

gas temperature, gas total pressure at the nozzle Inlet, the thrust reading, 
and tank pressure. These values, along with gas type, thrust calibration, the 
nozzle throat and exit diameter, and the number of data points were Input to a 
data reduction program. 

The data reduction program converted all the Input values Into SI units 
and calculated the pressure ratio, thrust, specific Impulse, thrust coeffi- 
cient, Reynolds number, discharge coefficient, and nozzle specific Impulse 
efficiency. The basic equations used In the data reduction are Included In 
appendix A. 


RESULTS AND DISCUSSION 

An evaluation of the performance differences of several nozzle contours 
was undertaken In this Investigation. The nozzles under consideration were for 
low thrust electrothermal devices which operate In the low Reynolds number 
range. In this flow regime, the major loss mechanism Is due to viscous effects 
which are associated with the boundary layer of the flow. Host of the past 
studies were performed with a heated gas flow, which made It difficult to eval- 
uate the loss mechanisms. Therefore, In an attempt to better understand these 
mechanisms, gases at ambient temperature were used. Thus, the heat transfer 
loss mechanism due to hot gas flow was eliminated. 

The experimental portion of this study evaluated the performance of noz- 
zles with conical, bell, and trumpet diverging contours. The parameters of 
Importance In the study of unheated flow Include the specific Impulse, the 
specific Impulse efficiency, the discharge coefficient, and the thrust coeffi- 
cient. The specific Impulse efficiency, specific Impulse, and thrust coeffi- 
cient are performance type parameters which gave similar results In nozzle 
performance trends; therefore, only the results of one of these parameters was 
Included. The experimental results Included were specific Impulse efficiency, 
and discharge coefficient variations with Reynolds number. 

The numerical analysis was used only to evaluate the conical nozzles 
because of the limitations of the numerical code. In essence, this code was 
run to determine where the slender-channel approximation was valid for the noz- 
zles considered and also to compare the predicted results with the experimen- 
tally measured performance. 
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Experimental Results 


Four variables were Involved In the nozzle testing; therefore, the 
approach to the data presentation had to take all four variables Into consider- 
ation. The data presented considers each gas separately over a range of 
Reynolds numbers. The Reynolds number range for nitrogen was from 500 to 6000 
and the range for hydrogen was from 150 to 3500. The data was plotted to 
determine the performance variation for each contour and for each area ratio. 
Several loss mechanisms were considered In the study, Including losses due to 
Incomplete expansion, viscous losses, and divergence losses. 

The specific Impulse (I sp ) efficiency Is defined as the ratio of actual 
specific Impulse to the theoretical maximum x sp; This efficiency Indicates 
the nozzle losses due to viscous and divergence losses. The theoretical maxi- 
mum specific Impulse of nitrogen and hydrogen at room temperature are 80 and 
300 sec, respectively. The data obtained for I sp efficiency was plotted 
for each each nozzle contour in order to observe the effect of area ratio. It 
was also plotted for each area ratio in order to compare the variation due to 
contour differences. Each r sp efficiency plot also shows the effect of 
Increasing Reynolds number. 

The l sp efficiency for different area ratio conical, bell, and trumpet 
nozzles run on nitrogen are shown In figures 9(a) to (e). As can be seen, the 
l sp efficiency decreases with decreasing Reynolds numbers for all the nozzles. 
The maximum difference In efficiency for the different area ratio 15° conical 
nozzles was 3 percent at the lowest Reynolds number and 2 percent over the rest 
of the range as shown In figure 9(a). The 20° conical nozzles, figure 9(b), 
had a 2 percent variation In the .Jsp.. efficiency over the entire Reynolds 
number range. The data for both then 0 and 20° conical nozzles were within 
expected experimental error; therefore, no area ratio could be considered sig- 
nificantly better than another. However, for the 25° conical nozzles, there 
does appear to be some difference In the performance for different area ratios. 
As the Reynolds number Increased, the efficiency of the nozzle with an area 
ratio of 25:1 did not Increase to the same level as the other 25° nozzles, 
figure 9(c). Its efficiency was about 3 percent lower than the other nozzles 
at the higher Reynolds numbers. The efficiency of the conical nozzles were all 
within 5 percent of each other over the entire Reynolds number range. Although 
differences were noted for these nozzles, the differences were within experi- 
mental error for the low Reynolds number cases, and just outside experimental 
error for the higher Reynolds number cases. The error bars on figure 9(a) 
Indicate this. 

The efficiencies of the bell contour nozzles ranged from 72 to 87 percent 
as shown In figure 9(d). At the lower Reynolds numbers, the nozzles with an 
area ratio of 50:1 had slightly higher efficiency. The data showed higher dif- 
ferences at the low Reynolds numbers and area ratio appeared to have little 
effect as the Reynolds numbers was Increased. 

For the trumpet nozzles, the efficiency differences were about 4 percent 
below a Reynolds number of 2000, and about 3 percent at higher Reynolds num- 
bers. The 25:1 area ratio nozzle had the lowest efficiency over the entire 
Reynolds number range. Again, at low and high Reynolds numbers the differences 
due to area ratio were within experimental error. 
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When the data for hydrogen was presented In the same fashion, the effi- 
ciency trends were similar for each group of nozzles. However, the effect of 
area ratio was more pronounced compared to the nitrogen data, as seen In 
figures 10(a) to (e). For the conical nozzles, with half-angles of 15°, 20°, 
and 25°, figures 10(a) to (c), respectively, as the Reynolds number Increased 
above about 1500, differences In efficiency with area ratio were apparent. The 
efficiency Increased with Increasing area ratio; the efficiencies of the 200:1 
area ratio nozzles were approximately 4 percent higher than those of the 25:1 
nozzles for both the 20° and 25° conical nozzles, and 3 percent higher than the 
15° conical nozzles with a 50:1 area ratio. The efficiency Increased about 
1 percent with each Increase In area ratio. Data presented by Murch (ref. 8) 
for 20° conical nozzles run with hydrogen at 815 °C (1500 °F) showed the oppo 
site trend of efficiency with area ratio. The percent differences In the Murch 
data were not discussed making It difficult to determine whether this 1 percent 
difference fell within experimental error. The efficiency was found to 
Increase with decreasing area ratios ranging from 20:1 to 200:1. The losses 
due to viscous flow can be Increased due to heat transfer at the nozzle wall, 
and these losses would be expected to Increase with higher area ratio and 
therefore, longer nozzles. However, because of the experimental error In the 
data presented here, and the unknown value of the experimental error In the 
Murch data, It Is possible that the data are consistent. 

The bell nozzles tested with hydrogen had efficiencies generally lower 
than the conical nozzles discussed above. However, the efficiency generally 
Increased with Reynolds number at lower area ratios (fig. 10(d)). The nozzle 
with the 25:1 area ratio had the highest efficiency over the range of Reynolds 
numbers tested. The differences In the efficiency with area ratio tended to 
decrease with Increasing Reynolds number. This trend was expected because 
these nozzles were designed assuming Isentroplc flow. The efficiency is 
expected to peak at higher Reynolds number than tested In this study, and the 
efficiency of the bell contoured nozzles at higher Reynolds number may be 
larger than any of the other nozzles considered here because as Reynolds number 
Increases the conditions can approach Isentroplc flow. 

The trumpet nozzles had efficiencies In the same range as those of the 
conical nozzles, with the efficiencies In the Reynolds number range below 2000 
slightly higher than the efficiencies of the conical nozzles, as shown In 
figure 10(e). The trumpet nozzle efficiencies only Increased slightly at 
Reynolds numbers above 2000. The larger area ratio nozzles were generally 
higher In efficiency over the Reynolds number range for hydrogen. 

Figures 11(a) to (e) show the efficiency differences of the contours at 
each area ratio for the tests run with nitrogen. For Reynolds numbers below 
2000, the bell nozzles had the lowest efficiency of all the area ratios tested. 
Figure 11(a) shows the efficiencies for the 25:1 area ratio nozzles. Below a 
Reynolds number of 2000, the efficiencies of the 25:1 conical and trumpet noz 
zles were within experimental error of each other. The efficiency of the bell 
nozzle was lower and dropped to 5 percent below the other contours at a 
Reynolds number of 500. The 15° and 20° conical nozzles had efficiencies 2 to 
4 percent higher at Reynolds numbers above 2000 for area ratios of 50:1 and 
25:1, respectively. For area ratios of 100:1 and greater, shown In 
figures 11(c) to (e), all the data was within 3 percent at Reynolds numbers 
above 2000. The only trend among these different higher area ratios was that 
the efficiency of the bell nozzles was lower for each. Also, the bell effi- 
ciency was still increasing while the other efficiencies remained constant. 
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This was expected because of the design technique used for the bell nozzles. 
That Is, at higher flowrates the flow Is more Isentroplc. At lower Reynolds 
numbers the trumpet nozzles generally gave higher efficiencies than those of 
the other contours. 

Data for the various nozzles operated on hydrogen are shown In 
figures 12(a) to (e). The data for area ratios of 25:1 and 50:1, presented In 
figures 12(a) and (b), display little difference between the performance of the 
different contours. At Reynolds numbers below about 1000, as the area ratio 
was Increased, figures 12(c) to (e), the 25° conical and trumpet nozzles per- 
formed slightly better than the other nozzles. There was likely some trade-off 
between the divergence losses In these more rapidly expanding nozzles and the 
viscous losses along the nozzle wall, as the area ratio was increased. 

The efficiency of the bell nozzles In the low Reynolds number range was 
lower than the other nozzles and dropped with Increasing area ratio. The bell 
nozzle designs were not optimized to Include boundary layer effects, as is done 
with larger thrust chemical propulsion nozzles where the boundary layer Is 
known to be thin. In these low Reynolds number nozzles, the flow was viscous 
with a thick boundary layer and, as found In studies by Rothe (ref. 1), may 
become fully viscous with no Indication of any Inviscld core for Reynolds num- 
ber less than 300. It was decided that optimization of the bell design for 
boundary layer effects was not practical for flows with thick boundary layers 
or for fully viscous flows. The comparison between different nozzle contours 
was on the basis of geometrically similar area ratios, which had differences 
of 4 percent or less at the same area ratio. 

Generally, the efficiencies of the nozzles using either nitrogen or hydro- 
gen were within experimental error of each other. The error was 5 percent for 
nitrogen data below Reynolds numbers of 2000 and for hydrogen data below 
Reynolds numbers of 1000. At higher Reynolds numbers, the error was 2 percent 
for both gases. The trends In both sets of data were similar. The bell noz- 
zles as designed were lower In efficiency, except for the 25:1 area ratio, 
which was within the scatter of the other data. Efficiency of bell nozzles 
would change If the boundary layer was taken Into account, as was done In an 
Investigation by Brophy, et al. (ref. 9), on arcjet nozzles. The trumpet noz- 
zle efficiency appeared to be better for Reynolds numbers below 2000, particu- 
larly at higher area ratios. However, given the range of the accuracy of the 
thrust stand, It appears that there was little difference between the nozzle 
configurations. The efficiencies of the 25° conical nozzles were generally 
within 1 to 2 percent of the efficiencies of the trumpet nozzles. This Is well 
within experimental error. It should be further noted that the conical nozzles 
have the advantage that they are typically simpler to fabricate. 

The data presented by Murch (ref. 8) for tests with hydrogen showed that 
as Reynolds number Increased, the bell and the 20° conical nozzles with area 
ratios of 100:1 were within 1 percent of each other, with the bell being the 
larger. The data for the trumpet nozzles showed that Its efficiency was about 
1 percent above that of the bell and conical shapes for Reynolds numbers of 
between 400 and 1200. Murch did not Include design details or detailed test 
procedures; therefore, the comparison can only be made on a relative basis. 

The trends exhibited were similar, although the efficiency values were 
different. 
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Recently Grlsnik, et. al. (ref. 11), reported results of a study In which 
unheated nitrogen and hydrogen were used to test a 20° conical, a bell, a 
trumpet, and a modified trumpet nozzle. This data showed all of the nozzles 
tested had the same performance to within experimental error, which was 
5 percent In that study. It should be noted that nozzles compared by Grlsnik 
had area ratios ranging from 120:1 and to 150:1. The nozzle efficiency ranges 
obtained for the hydrogen data were similar to the data obtained In this study. 
The nitrogen data reported was 3 to 4 percent higher than the data obtained In 
this study at Reynolds numbers below 2000, but the trends shown were similar. 

The discharge coefficient was the other parameter selected for study In 
this experimental Investigation. The discharge coefficient Is defined as the 
ratio of the measured mass flowrate to the theoretical maximum mass flowrate 
through a nozzle assuming Isentroplc flow. The discharge coefficient Indicated 
the mass losses due to viscous effects In the nozzle throat. The discharge 
coefficient data was presented In the same format as the specific Impulse effi- 
ciency data with variation of contour and area ratio presented as a function 

of Reynolds number. 

The nitrogen discharge coefficient data for the different nozzle contours 
are presented In figures 13(a) to (e). Results for the 15°, 20°, and 25° coni- 
cal nozzles are shown In figures 13(a) to (c). As can be seen, the flow coef- 
ficient decreased with decreasing Reynolds number. This trend was expected 
because at lower Reynolds numbers the boundary layer extends across more of the 

throat and, thus, more of the mass flows through the viscous layer. As the 

Reynolds number Increased, the discharge coefficient rose to a constant value 
and showed no difference between area ratios, which was also expected. 

The discharge coefficient data for the five bell contoured nozzles Is pre- 
sented In figure 13(d). The discharge coefficients are again seen to be the 
same within the error of the flow measurement at the higher Reynolds numbers. 

At low Reynolds numbers, the discharge coefficient dropped as expected, except 
for the nozzle with the 25:1 area ratio. It Is suspected that differences In 
the throat curvature may have accounted for this performance difference, but 
the bell nozzles all have the same Inlet throat radius of curvature and the 
same, although different from the inlet, outlet curvature. Another possibility 
for discharge coefficient not behaving as expected may be surface roughness 
effects In the throat region. However, surface roughness was expected to have 
a negligible effect with the thick boundary layers occurring at these low 
Reynolds numbers. 

The trumpet nozzle data for discharge coefficient Is presented In 
figure 13(e). The trumpet nozzle discharge coefficients were lower than the 
other contours, at low Reynolds numbers and the percent difference was about 
4 percent at the higher Reynolds numbers. The lower discharge coefficient was 
an Indication that the throat boundary layer In the trumpet nozzles was larger, 
thus, the viscous losses In the throat region were higher. 

The discharge coefficient data for the different nozzle contours run with 
hydrogen Is presented In figures 14(a) to (e). The data In figures 14(a) 
to (c) showed the same trends as the discharge coefficient data of the 15°, 

20°, and 25° conical nozzles tested with nitrogen. However, as the Reynolds 
numbers Increased, the discharge coefficient values did not Increase to the 
same level as those for the nitrogen data. For Reynolds numbers higher than 
600, the flowrate data was obtained using the flowmeter with a range of 0 to 
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50 slpm described In the Experimental Apparatus section. This flowmeter was 
found to consistently read 3 to 5 percent low over the range used (6 to 
26 slpm). Since the discharge coefficient Is a ratio Involving the actual mass 
flowrate, It follows that the values obtained would likewise be expected to be 
3 to 5 percent low. Correction of the data by 3 to 5 percent for Reynolds 
numbers from 600 to 3600 would raise the level of values up to the same range 
as obtained In the nitrogen tests. The nitrogen data spanned a flowrate range 
of 0.5 to 6 slpm, and used a flowmeter with an error of less than 1 percent. 

Figure 14(d) shows the discharge coefficients of the bell contour nozzles. 
Again, the data followed the expected pattern except for the 25:1 area ratio 
nozzle. It was not clear why this discrepancy occurred. The data above a 
Reynolds number of 600 should be corrected upward by 3 to 5 percent as a result 
of error In the flow measurement. 

At Increasing Reynolds numbers the discharge coefficients are less dis- 
tributed and appear to approach some limiting value. The trumpet nozzle dis- 
charge coefficient data for hydrogen, shown In figure 14(e), was lower over the 
entire range of Reynolds numbers than the other nozzles. Correction of the 
measured flowrate would, again, raise the data for Reynolds numbers above 600 
to the values obtained with nitrogen. The Implication again was that for the 
lower Reynolds numbers the viscous losses In the throat region were greater 
than those of the other contours. 

Figures 15(a) to 15(e) are plots of the discharge coefficients for nitro- 
gen for different nozzle contours at each area ratio. At lower area ratios of 
25:1 and 50:1, shown In figures 15(a) and (b), respectively, the trumpet noz- 
zles had lower discharge coefficients. As the area ratio Increased, the 
trumpet nozzles had lower discharge coefficients In the low Reynolds number 
range, as shown In figures 15(c) to (e) for area ratios of 100:1, 150:1, and 
200:1, respectively. The discharge coefficients of the bell nozzles also 
tended to decrease with Increasing area ratio for the lower Reynolds numbers. 

At larger Reynolds numbers, the 15°, 20°, and 25° conical nozzles had essen- 
tially the same discharge coefficients. When this data was plotted for hydro 
gen, the same general trends were exhibited, as shown In figures 16(a) to (e) 
for area ratios of 25:1, 50:1, 100:1, 150:1, and 200:1, respectively. The 
variations In discharge coefficient were attributed to differences In boundary 
layer thickness In the throat region, and viscous losses associated with the 
boundary layer. 

Several Ideas were explored as a means of explaining the variation In dis- 
charge coefficient results, among the possible factors was the Influence of 
geometry In the throat region. The geometry of the throat region Included the 
converging section, the throat radius of curvature, and the surface condition 
of the nozzle wall. As mentioned earlier, the Inlet side or convergent section 
of each of the various nozzle contours was the same. The convergent angle was 
45° In each nozzle and the measured throat diameters were all within 2 percent 
of each other. Accordingly, the Inlet geometry of the different nozzles was 
not considered to be a factor In the variation of the discharge coefficients. 

In addition, previous experimental and numerical studies, have Indicated that 
the convergent nozzle contour does not affect the nozzle discharge coefficient 
(refs. 20 and 21). 

Nozzle throat geometry begins to affect the flowfleld at the point of tan- 
gency between the convergent angle and the circular radius at the throat. For 
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the nozzles tested, the radius of curvature was designed to be the same on the 
throat Inlet side; however, the throat radius of curvature on the outlet side 
was different for each contour. For comparison each outlet radius of curvature 
was normalized to the throat radius. These ratios were 8 for the 15° conical 
nozzles, 6.7 for the 20° conical nozzles, 5.2 for the 25° conical nozzles, 1 
for the bell nozzles, and 40 for the trumpet nozzles. The ratios for the con- 
ical nozzles had to be estimated because the actual radius was obtained by 
blending the throat and divergent section contours. This was done to simplify 
the fabrication of these nozzles, hence the actual value could not be deter- 
mined. The trend In both the nitrogen and hydrogen tests was for Increases In 
discharge coefficient to correspond to decreases In the outlet radius of curva- 
ture. Data presented by Back (ref. 18) Indicated an opposite trend, l.e., dis- 
charge coefficient Increased with Increasing radius of curvature. However, two 
significant differences exist between these sets of data. The data presented 
here was obtained with nozzles that did not have a constant radius of curvature 
through the throat, and the data was obtained for low Reynolds number flows 
(<10 000). On the other hand, Back's data was for high Reynolds number flow 
through nozzles with a constant radius of curvature through the throat. These 
differences make It difficult to make valid comparisons, but they Indicate that 
this may be an area for further study. An Important point here was that 
despite the low discharge coefficients In the low Reynolds number region, the 
trumpet nozzles, at least at higher area ratios, appeared to have better I sp 
efficiencies than the other nozzles. Apparently, the gains made In the diver- 
gent section of the trumpet nozzles overshadowed the viscous losses In the 
throat region. It was suspected that further study of the throat region may 
further Improve the performance of these small nozzles. Through Improved 
design of the throat, the trumpet design may obtain even better efficiencies. 

The surface condition of the nozzle throat region after machining may 
possibly affect the discharge coefficient. However, this was not considered 
likely to be significant for the lower Reynolds number flows. In low Reynolds 
number flows, the boundary layer was thick, as Indicated by Rothe (ref. 1) and 
others. Thus, the surface roughness of the wall In the throat region and down- 
stream were not likely to have a significant effect on the flow field unless 
the surface roughness was on the order of the boundary layer thickness, as Is 
possible In high Reynolds number flow. The effect of surface roughness should 
be verified In subsequent testing, but, as Indicated above. Is not expected to 
be significant In the low Reynolds number regime. 


Numerical Results 

A computer code for calculating viscous flow In nozzles was developed by 
Rae (ref. 25). The performance predictions given by this code were used for 
comparison with the experimentally measured performance. The method employed 
was based on the slender-channel approximation of the Navler-Stokes equations. 
The slender-channel equations are formally Identical to the boundary layer 
equations. However, unlike boundary layer solution procedures, the axial pres- 
sure gradient Is unknown and must be calculated as part of the solution. The 
equations are written in cylindrical coordinates for steady, axlsymmetrlc flow. 
The basic equations and their nondlmenslonal form are listed In appendix B. 

In the solution of the equations, the boundary conditions on the axis resulted 
in all radial derivatives vanishing. At the nozzle wall, rarefied flow bound- 
ary conditions were used, allowing a velocity slip and a temperature jump to 
exist. In this approach, Rae assumed the gas to obey the perfect-gas law, with 
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a viscosity that Is proportional to a power of temperature. The Prandtl number 
and the specific-heat ratio were taken to be constant. 

The numerical procedure began at a location upstream of the nozzle throat. 
For the cases considered here, this location was set at the beginning of the 
convergent section. The slender-channel equations were solved using an 
Implicit finite difference scheme, which marched downstream through the nozzle. 
The required Inputs to this code were: the reservoir conditions, nozzle geom 

etry, gas properties, and, If heat transfer was Included, the distribution of 
wall temperature. For the cases calculated In this study, the adiabatic wall 
boundary condition was used. This was considered a good approximation because 
the experimental tests selected for comparison were run with an unheated gas. 

This program has an option that permits the user to input either a known 
mass flowrate or Reynolds number wherein the flowrate Is calculated for the 
given nozzle geometry. All the cases run with nitrogen had the mass flowrate 
calculated as part of the solution. When tried with the known mass flowrate, 
some of the hydrogen and all of the nitrogen cases failed to converge to a 
solution. Therefore, for the data presented here only the nitrogen cases with 
calculated flowrates were Included. The hydrogen cases Included some with a 
known mass flowrate and some with code calculated mass flowrates. 

The Rae code was written for nozzles with conical geometries; therefore, 
only conical nozzle performance could be compared. For this comparative study, 
the 15°, 20°, and 25° half angle conical nozzles with area ratios of 50 and 150 
were selected. Data from these six nozzles at three different Reynolds numbers 
for both nitrogen and hydrogen were used as Input to the code. In order to run 
the program, the nozzle geometry had to be defined with axial length normalized 
to the throat radius. The geometry Included the half angle of the converging 
and diverging sections, the throat radius of curvature and definition of the 
axial length. These values are listed In table I for each nozzle considered. 
The throat of each nozzle was defined as the axial origin, the convergent sec- 
tion was then In the negative axial direction and the divergent section 
extended In the positive direction. The code was built on the assumption that 
the throat radius of curvature was constant; however, the conical nozzles were 
designed such that the radius of curvature on the Inlet side of the throat was 
different from that on the outlet side. The Inlet side of each nozzle had the 
same design; therefore, the Inlet curvature was used. A complete list of Input 
variables Including their definitions Is given In appendix C. 

The numerical code developed by Rae (ref. 25) was used to generate per- 
formance predictions which were compared to the experimental performance data. 
The results obtained from the code for nitrogen and hydrogen are listed In 
tables II and III, respectively. The tables list the nozzles area ratio and 
half angle, the Reynolds number (B), the experimental flowrate (A code), the 
flowrate at which the code ran (Ag^ ven ), the centerline Mach numbers at the 
throat and exit, and the thrust coefficient. For each nozzle and gas, three 
different Reynolds numbers flows were examined. In both the nitrogen and 
hydrogen cases, the highest Reynolds number case treated did not converge to a 
solution for any of the nozzles. For nitrogen, the Reynolds numbers were 
approximately 900, 3500, and 8900. The Reynolds numbers for hydrogen were 
about 550, 3750, and 6800. These Reynolds numbers have been redefined as 
required prior to Input to the code. 
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The nitrogen cases listed In table II only ran to a solution when the code 
calculated the flowrate (A), and resulting flowrates were, with the exception 
of run number 20, lower than the experimentally measured values. The flowrates 
that were calculated are proportional to discharge coefficient and are related 
to Reynolds number (B). The Reynolds number Input to the program was calcu- 
lated as a function of nozzle Inlet stagnation conditions (enthalpy, tempera- 
ture, and pressure), and was not calculated the same way as the experimental 
Reynolds number. The latter was determined from the experimental mass flowrate 
and the throat conditions. Therefore, there was some mismatch or built-in 
error between the code required flowrate and Reynolds number. These values 
were not directly related as the experimental flowrate and Reynolds number 
were. This was believed to be the source of difficulty In running some of the 
cases with known flowrates. The code provided the option of calculating or 
Inputting a known flowrate. The cases with N2 were all code generated 
flowrates . 

For the hydrogen runs listed In table 111, several cases were found to 
converge to solutions with the given flowrates. Apparently, an error was 
Introduced Into the Reynolds numbers and the flowrates for the other cases, 
similar to that of the nitrogen data. Therefore, these cases were run allowing 
the code to calculate the flowrate. Several cases for both gases were unable 
to converge to solutions with either the known or code calculated flowrate. 

In Williams discussion of the slender-channel approximation (ref. 26), the 
transport properties were assumed to vary as a power of temperature. Williams 
Indicated that for a given wall angle and wall temperature, solutions were pos- 
sible only for certain conditions of throat centerline Mach number, Reynolds 
number, ratio of specific heat and Prandtl number. Rae (ref. 2) also Indicated 
that his code, which used this same approximation, had limitations and that 
there were concerns about accuracy. At low Reynolds numbers, It was probable 
that no solution with a supersonic core would be found. At high Reynolds num 
bers, accuracy could be reduced for two possible reasons: first, as the bound- 

ary layer thins, It would be described by fewer radial grid points, and second, 
the slip velocity decreases at the wall; thus, the Integrals that are related 
to velocity would not be calculated correctly by Simpson's rule. 

The cases that ran using the slender-channel approximation were compared 
to the experimental data. The performance parameter obtained from the code was 
thrust coefficient. The thrust coefficient Indicates the Increase In thrust 
due to expanding the gas through a nozzle as compared to thrust obtained If the 
chamber pressure acted only over the throat area. 

Ihe results obtained for the nitrogen cases are shown In table IV. The 
table gives the nozzle, the Reynolds number, the experimental and numerical 
thrust coefficient, and the percent difference between these values. The 
numerical results are for Reynolds numbers of 460 and 1850 (B of 900 and 3500), 
at area ratios of 50:1 and 150:1. The results showed that for 15° conical noz- 
zle, the low Reynolds number point for the 50:1 nozzle was 7.6 percent higher 
than the experimental value, and the other numerical thrust coefficients were 
3 to 3.5 percent higher. The results for the 20° conical nozzles showed the 
numerical results to be within 4 percent of the experimental thrust coeffl 
clents. The 25° conical nozzle results were within 3 percent of the experimen- 
tally obtained thrust coefficients. 
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The results obtained for hydrogen are shown in table V. Numerical results 
were obtained for Reynolds numbers of 260 and 1850 (B of 550 and 3750). In the 
hydrogen data, all the cases at a Reynolds number of 260 and the cases for the 
15° cone with an area ratio of 150:1 and a Reynolds number of 1850 had code 
generated flowrates. In these cases, the percent difference between the exper- 
imentally and numerically obtained thrust coefficients ranged from 5 to 
20 percent. Cases run with known flowrates showed differences of 3 percent or 
less. Kallis, et al. (ref. 29), used Rae's code to predict resistojet perform- 
ance and also compared the predictions to experimental data. The code was run 
for a range of Reynolds numbers from 1000 to 7000. Differences between experi- 
mental and predicted performance ranged from 0 to 9.2 percent. The predictions 
obtained here were similar, except for the hydrogen cases with code generated 
flowrates. 

The agreement between experimental thrust coefficients and computed thrust 
coefficients obtained from code generated flowrates for the two gases was very 
different. The nitrogen results agreed well, but the hydrogen agreement was 
poor in some cases. The poor results obtained for hydrogen were only for the 
low Reynolds number cases. It is suspected that the code was not accurate at 
modeling the Reynolds number range where most or all of the divergent section 
of the flow field was subsonic. As indicated earlier, Rae (ref. 25) had not 
accurately defined the limits of this code. 

After running the code, it was found that there appeared to be a numerical 
stability problem related to the throat centerline Mach numbers. It was found 
that the program would terminate at centerline Mach numbers of 1.2 or higher. 
Other cases were unable to calculate the pressure gradient within the required 
number of iterations, and some cases could not get beyond the throat. The 
cases that did not run stopped for one of these reasons. It Is believed that 
the higher Reynolds number cases that did not run could also have had geometry, 
Reynolds number or gas property conditions for which the slender-channel 
approximation was inappropriate. 

An interesting result in the hydrogen numerical results was that the 
centerline Mach number at both the throat and exit were generally higher for 
the low Reynolds number case. The velocity profiles at several locations along 
the nozzle from the throat to the exit are compared in figures 17(a) and (b). 
The two cases presented are for the 15° conical nozzle with a 50:1 area ratio, 
at Reynolds numbers of 260 and 1850. At the throat, it is seen that the veloc 
1 ty profiles are approximately the same. But, at axial locations farther down 
stream, the velocity close to the centerline (ETA = 0 to 0.2) was significantly 
higher for the low Reynolds number case. This suggests that the higher center- 
line Mach number, especially at the exit, was to be expected. 

Although the results predicted by Rae's code were reasonable in many 
cases, a word of caution is in order. As previously mentioned, the geometry 
written into this code and used in the cases presented here did not accurately 
describe the nozzles that were experimentally tested. Rae's code assumed a 
constant throat radius of curvature; however, the nozzles tested had different 
curvatures on the inlet and outlet sides of the throat. It is recommended that 
the geometric section of the code be refined to accommodate the variation of 
radius of curvature and to accommodate different contours. The addition of the 
throat curvature variation is also recommended because as previously discussed, 
some experimental studies have shown some differences in the discharge coeffi- 
cient with changes in the throat radius of curvature. 
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CONCLUSIONS 


The major goal of the study presented here was to evaluate and compare the 
performance of different nozzle contours. The study evaluated nozzles with 
different divergent section contours, Including 15°, 20°, and 25° conical noz- 
zles, bell nozzles, and trumpet nozzles. Area ratios of 25:1, 50:1, 100:1, 
150:1, and 200:1 were tested for each nozzle, except for the 15° conical nozzle 
In which all of these area ratios were tested except the 25:1. The nozzles 
were tested with unheated nitrogen and hydrogen over Reynolds number ranges of 
500 to 6000 and 150 to 3500, respectively. 

The performance of the nozzles was generally within experimental error of 
each other over the Reynolds number range tested. The trends seen In the spe- 
cific Impulse efficiency were similar for both gases. The bell contour nozzles 
tended to be lower In efficiency than the other contours In the low Reynolds 
number range tested. At Reynolds numbers below about 2000, the trumpet and 25° 
conical nozzles performed slightly better than the other conical nozzles, and 
significantly better than the bell nozzles, for higher area ratios. But, these 
results were still within experimental error, which makes any recommendation 
quite difficult. The discharge coefficient results showed that the trumpet 
nozzles had the lowest values, particularly at the lower Reynolds numbers. 

This may Indicate that the better performance of the trumpet nozzles may be 
further Improved through a more careful design of the throat region. Reduction 
of the viscous losses at the throat could Increase the overall performance of a 
nozzle. 

Although the trumpet or 25° conical nozzles appeared to be slightly better 
designs, at lower Reynolds numbers, It Is unclear which nozzle Is better as all 
fell within experimental error. It Is recommended that a similar comparison be 
made between these nozzles tested with hot flow to evaluate the effects of heat 
transfer on overall performance. The bell nozzle as designed Is not recom- 
mended for use In low Reynolds number nozzles. Consideration should also be 
given to fabrication of these small nozzles. If performance differences are 
small, as Indicated by these tests, and not a major driver In a particular 
application, the choice may be based Instead, on which nozzle Is easiest to 
fabricate. 

The numerical results obtained using a computer code based on the 
slender-channel approximation, showed reasonable agreement with experimental 
results for nitrogen and for some of the hydrogen data. However, the code used 
did not correctly define the geometry of the nozzles tested. The code should 
be modified to more correctly define the nozzle geometry and therefore Improve 
confidence In the performance prediction. The code was also found to have 
limited use In the Reynolds number range considered here, as no solutions were 
obtainable for the highest Reynolds numbers tried. 
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APPENDIX A - BASIC EQUATIONS FOR DATA REDUCTION 


The measured thrust based on a momentum balance at the nozzle is 

F m = roU E + A E (P E " V 

For the same flow conditions, the thrust in hard vacuum is 

Fy = mU E + P e A e (/ 

To account for the effect of the vacuum tank pressure in excess of hard 
vacuum, thrust is 

F. = F + P.A (/ 

A m A E 

Specific impulse was then obtained using F/\; 


(rti g) 

Thrust coefficient determined the Increase in the thrust due to expanding 
the gas through a nozzle compared to thrust if the chamber pressure acted over 
the throat area only. It was calculated by; 


T " (P C A*) ' ' 

The Reynolds number was calculated based on throat conditions and was 

_ 7-^rr (A6) 

R = (ify d) 

The discharge coefficient indicated the mass flow losses due to viscous effects 
in the nozzle throat. It is the following ratio: 


The nozzle specific Impulse efficiency is 


i _ T sp,A 
sp sp,I 
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LIST OF SYMBOLS 


Af; nozzle exit area, m 2 

A* nozzle throat area, m 2 

Cq discharge coefficient 

Cj thrust coefficient 

d throat diameter, m 

F a adjusted thrust, N 

F m measured thrust, N 

g 9.807 m/sec 2 

I S p specific Impulse, sec 

I S p,A measured specific Impulse, sec 

I sp * I specific Impulse assuming Isentroplc flow, sec 

m mass flowrate, kg/sec 

it)A measured mass flowrate, kg/sec 

mj theoretical maximum mass flowrate, kg/sec 

N. specific Impulse efficiency 

sp 

Pa ambient (tank) pressure, N/m 2 

P c Inlet pressure, N/m 2 

P£ static pressure at exit plane of nozzle, N/m 2 
R Reynolds number 

U E average nozzle exit velocity, m/sec 

v viscosity, kg/(m sec) 
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APPENDIX B - BASIC EQUATIONS OF RAE'S CODE 

For flow In slender-channels , the approximate equations of motion for 
steady axlsymmetrlc flow In cylindrical coordinates are: 

continuity, 

fj (pu) + (pv) + Jr = 0 (Bl) 

axial momentum, 

9u au dp la/ au\ 

M * ' ,v ■ - dZ * Tr) < e2 > 

radial momentum, 



( B3) 


energy, 


ah 


ah 


iz + pv if 



(B4) 


The coordinates and the dependent variables are made dimensionless using 
throat radius and the reservoir conditions as: 


x 



U 


u 


a = 


r 

r* 


V 


■ 


p » 




In terms of these variables, the equations of motion become: 
continuity. 


momentum, 



PnU _ 
we* 


£ /„ au w_ au\ y - l d£ _l_ i_ ( au\ 
0 \ 3x % *") 2y dx Bno 2 an \ n an ) 


(B5) 


(B6) 
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energy. 


£ l ii £© . w_ m\ . y - 1 dP J a_ w ae\ 2e^ ( au\ 2 

0 V 3X °W 3n / = Y dx Bno 2 3n V r 3 V Bo* W 


(B7) 


In these equations, density was eliminated by using the equation of state: 
P = Oe. The radial velocity was transformed and appears at W = V - Un do u /dx. 


LIST OF SYMBOLS 

B Reynolds number based on reservoir conditions 

0 p/p 0 

H total enthalpy 

h static enthalpy 

P P/P 0 

p pressure 

Pr Prandtl number 

r,z cylindrical coordinates 

U.V u/V2i£, v/\/2U^ 

u,v axial and radial velocity components 

W V - Un do w /dx 

x,o Z/r*, r/r* 

y specific heat ratio 

n o/<? u 

8 h/H 0 

w viscosity 

p density 

o w r w /r* 

o> exponent In viscosity and enthalpy relation 

subscripts 

o at reservoir conditions 
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w at wall 

* at geometric throat 
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APPENDIX C - INPUT VARIABLES FOR VISCOUS FLOW CODE 

AO Initial lower boundary of flowrate 

A1 Initial upper boundary of flowrate 

ALPHT temperature accommodation coefficient 

ALPHU velocity accommodation coefficient 
ATES1 tolerance allowed for A 

B reservoir Reynolds number 

DELX Initial delta X 

GAMMA specific heat ratio 

OMEGA temperature exponent In viscosity law 
PR Prandtl number 

R1 throat radius of curvature normalized to throat radius 

THE1A1 half angle of convergent cone 

THETA2 half angle of divergent cone 

TW2 wall to temperature ratio (heat transfer case) 

XO Initial X value (negative) 

XI begin reduced step size (throat region) 

X2 restore original step size 

XCUT switch point for straight line extrapolation through saddle point 

XIPG defines method of calculating pressure gradient 

XMAX calculations stopped at this point 

XPRIN1 profiles printed at this Interval 

XTW1 used for heat transfer case 

XTW2 used for heat transfer case 

IBC =l,heat transfer; =2, adiabatic wall 

1ETAPR Interval at which profiles printed 

ISP -0 , print Iterations; =1, print profiles only 

MORE If nonzero, more than one case will be run 
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NPRO 

NRUN 


If nonzero, profiles are not printed 
run number 
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TABLE I. - DEFINITION OF NOZZLE 
GEOMETRY 



e T 

€ 

e 2 

r/r* 

o 

X 

X max 

4 

5 

50 
T 50 
50 
T 50 
50 
T 50 

T 5 
T 5 
20 
20 
25 
25 

2. 

67 

r 

-5 

TT .5 
20.9 
8.6 
T 5 . 5 
6.8 
12.2 


TABLE II. - NUMERICAL METHOD RESULTS FOR NITROGEN 


Run 

number 

C 

e 2 

B 

A gi ven 

A code 

M* 

M e 

C T 

20 

50 

15 

999 

tm 

0.1160 

0.962 

3.42 


21 

50 



3600 


.1239 

1 .002 

4.17 

1.524 

22 

50 



8823 

■EH 

- - - • ■ 

— 

— 


23 

150 



923 

.1257 

.1168 

1 .026 

4.65 

1.294 

24 

150 



3534 

.1287 

.1241 

1 .020 

5.07 

1.537 

25 

150 



8914 

.1281 

~ 


. 

... .. 

26 

50 


890 

.1277 

.1159 

.994 

3.74 

1 .364 

27 

50 



3536 

.1286 

.1241 

1.020 

4.05 

1.534 

28 

50 



8880 

.1280 


— . 


. 

29 

150 



944 

.1204 

.1159 

.976 

4.30 

1.378 

30 

150 



3569 

.1280 

.1241 

1 .017 

4.90 

1.556 

31 

150 



9029 

.1267 

— 

— 

— 



32 

50 

25 

918 

.1238 

.1159 

.984 

3.75 

1 .386 

33 

50 



3542 

.1283 

.1241 

1 .019 

4.06 

1.543 

34 

50 



8949 

.1278 


— 




35 

150 



923 

.1232 

.1168 

1 .026 

4.59 

1 .400 

36 

150 


1 

3560 

.1277 

.1241 

1 .018 

4.96 

1.569 

37 

150 


r 

8941 

.1279 

— 

— 


— 
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TABLE III. - NUMERICAL METHOO RESULTS FOR HYDROGEN 


Run 

number 

C 

e 2 

B 

Agl ven 

A code 

M* 

Me 

C T 

T 

El 

T 5 

561 

0.1154 

0.1120 


3.91 

1.238 


B9 



3746 

.1226 

.1226 

.932 

3.20 

1 .469 


50 



6556 

.1246 



— 



T 50 



531 

.1232 

.1168 

1 .192 

1 .73 

1.132 


T 50 



3760 

.1212 

.1241 

1.002 

4.99 

1.540 


1 50 



6850 

.1234 



— 

— 


50 


548 

.1188 

.1168 

1.178 

3.81 

1 .328 


50 



3724 

.1225 

.1225 

.929 

2.54 

1 .408 


50 



6801 

.1245 



— 

— 

TO 

T 50 



587 

.1109 

.1109 

.933 

4.21 

1.274 

TT 

T 50 



3788 

.1203 

.1203 

.865 


1 .496 

T 2 

T 50 



6919 

.1223 


— 



— 

T 3 

50 

25 

571 

.1141 

.1120 

.977 

3.65 

1.307 

T 4 

50 



3761 

.1215 

.1215 

.897 

3.33 

1 .424 

1 5 

50 



6853 



— 

— 


T 6 

T 50 



543 

.1199 

.1168 

1 .182 

4.41 

1 .338 

17 

T 50 



3790 


.1204 

.868 

3.92 

1 .496 

18 

T 50 



6896 

.1224 


— 

— 

MB 


TABLE IV. - COMPARISON OF NITROGEN 
EXPERIMENTAL AND NUMERICAL 
PERFORMANCE RESULTS 


e 2 

C 

Re 

m 

C T 

numer- 

ical 

Percent 

difference 

15 

50 

460 

1.29 

1 .341 

+7.6 

15 

50 

1845 

1.25 

1.524 

+3.1 

15 

150 

460 

1.25 

1.294 

+3.5 

15 

150 

184/ 

1 .50 

1.537 

+3.4 

20 

50 

458 

1.40 

1.365 

-2.2 

20 

50 

1830 

1 .51 

1.534 

+1.6 

20 

150 

460 

1.33 

1 .378 

+3.9 

20 

150 

1848 

1.52 

1.556 

+2.4 

25 

50 

458 

1.34 

1.384 

+2.9 

25 

50 

1848 

1.50 

1.556 

+3.0 

25 

150 

462 

1.38 

1.399 

+1 .3 

25 

150 

1848 

1.53 

1.569 

+2.5 
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TABLE V. - COMPARISON OF HYDROGEN 
EXPERIMENTAL AND NUMERICAL 
PERFORMANCE RESULTS 


e 2 



C 


m 

C T 

numer- 

ical 

Percent 

difference 

15 

50 

263 

mm 


+ 7.7 

15 

50 

1836 

wttm 

1.469 


15 

150 

265 

wRM 

1 .132 

-8.0 

15 

150 

1836 

1.47 

1.540 

+4.8 

20 

50 

265 

1.11 

1 .328 

+20 

20 

50 

1852 

1 .46 

1 .408 

+3.3 

20 

150 

262 

1.11 

1.276 

+14 

20 

150 

1838 

1 .46 

1 .496 


25 

50 

260 

1 .16 

1 .307 

+12 

25 

50 

1834 

1 .46 

1.474 

+1.1 

25 

150 

261 

1.26 

1.338 

+5.9 

25 

150 

1829 

1 .46 

1.496 

+2.3 
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FIGURE 2. - NOZZLE ASSEMBLY ON THRUST STAND. 
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C-87-2^75 

FIGURE 3. - DIFFERENT AREA RATIO 20° CONICAL NOZZLES WITH FLANGE. 


TABLE OF DIMENSIONS/ MM 

PART 

e 

IDE 

* 

£ 

L i 

L 2 

NUMBER 






6152 

15° 

10.79 

50 

17.27 

25.91 

6153 



15.26 

100 

25.60 

34.24 

6154 



18.64 

150 

31.98 

40.64 

6155 



21.53 

200 

37.36 

46.00 

6201 

20° 

7.66 

25 

8.38 

17.02 

6202 



10.77 

50 

12.70 

21.34 

6203 



15.14 

99 

18.85 

27.48 

6204 



18.65 

150 

23.55 

32.18 

6205 



21.56 

200 

27.50 

36.14 

6251 

25° 

7.69 

25 

6.53 

15.16 

6252 



10.78 

50 

9.93 

18.54 

6253 



15.14 

99 

14.70 

23.34 

6254 



18.67 

150 

18.39 

27.03 

6255 



21.50 

199 

21.49 

30.12 


*£ = AREA RATIO. 

THROAT DETAIL 


1.524 



APPROXIMATE R ^ 



FIGURE 4. - CONICAL NOZZLE DESIGN. (ALL DIMENSIONS ARE IN MILLI- 
METERS.) 
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FIGURE 5. - DIFFERENT AREA RATIO BELL NOZZLES WITH FLANGES. 


TABLE OF DIMENSIONS/ mm 


PART 

NUMBER 

IDE 

• 

e 

L i 

L 2 

8.61 

7.645 

26 

10.16 

18.80 

8.62 

10.87 

52 

16.40 

25.04 

8.63 

15.24 

104 

26.34 

34.96 

8.64 

18.62 

154 

34.70 

43.33 

8.65 

21.54 

205 

42.16 

50.80 


*£ = AREA RATIO. THROAT DETAIL 

0.762 



FIGURE 6. - BELL NOZZLE DESIGN. (ALL DIMENSIONS ARE IN MILLIMETERS.) 
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FIGURE 7. - DIFFERENT AREA RATIO TRUMPET NOZZLES WITH FLANGES. 


TABLE OF DIMENSIONS, mm 


PART 

NUMBER 

DIA 

L i 

L 2 

£* 

H425 

7.643 

9.14 

17.02 

26 

H450 

10.78 

10.95 

18.82 

51 

H4100 

15.24 

12.73 

20.60 

103 

mm 

18.66 

13.69 

21.56 

155 

m 2 oo 

21.58 

14.32 

22.20 

207 


THROAT DETAIL 


£ = AREA RATIO. 


1.83 R 



FIGURE 8. - TRUMPET NOZZLE DESIGN. CALL DIMENSIONS ARE IN MILLI- 
METERS.) 
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(B) CONICAL NOZZLES WITH 20° HALF ANGLE. 
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(D) BELL NOZZLES. (E) TRUMPET NOZZLES. 

FIGURE 9. - NOZZLE EFFICIENY VARIATION WITH AREA RATIO, UNHEATED NITROGEN. 
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FIGURE 15. - DISCHARGE COEFFICIENT FOR VARIOUS CONTOURS WITH UNHEATED NITROGEN. 
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FIGURE 17. - VELOCITY PROFILES AT VARIOUS AXIAL LOCATIONS FOR THE 15° CONICAL WITH AREA RATIO 50:1 
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